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Abstract 

qh| We study the infrared effective theory of gravity that stems from the quantum 

trace anomaly. Quantum fluctuations of the metric induce running of the cosmological 
bX), constant and the Newton constant at cosmological scales. By imposing the generalized 

Bianchi identity we obtain a prediction for the scale dependence of the dark matter and 
' dark energy densities in terms of the parameters of the underlying conformal theory. 

For certain values of the model parameters the dark energy equation of state and the 



I observed spectral index of the primordial density fluctuations can be simultaneously 

OO • reproduced. 

m ■ 



^ ■ 1 Introduction 



The evidence of the cosmic acceleration presumably driven by dark energy (DE) with neg- 
ative pressure [1] and precise measurements of the cosmic microwave background radiation 
[21 [3] have triggered a renewed interest in the cosmological constant Among different 
approaches to study the cosmological constant a convenient framework is based on the effec- 
tive field theory [5l[6]. This theory is a long distance realization of quantum gravity reduced 
to the general theory of relativity supplemented by the quantum field theory in curved space 

The cosmological constant A and the Newton constant G, being parameters in the 
Einstein-Hilbert action, receive contributions from quantum loops and become running con- 
stants, i.e., functions of the running scale [9l [TOl [11]. Unfortunately, the running scale, 
intuitively expected to be of the order of typical momenta of the particles in loops, cannot 
be fixed unambiguously. 

It has been shown [12j that a consistent infrared modification of gravity due to the 
quantum trace anomaly implies the presence of additional terms in the low energy effective 
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action. Quantum fluctuations of tlie metric modify the classical metric description of general 
relativity at cosmological scales and provide a mechanism for a screening of the cosmological 
constant and the inverse Newton constant. 

The effective theory of the conformal factor induced by the quantum trace anomaly has 
a non-trivial infra-red (IR) dynamics owing to the existence of a non-trivial IR stable fixed 
point [13j. A number of issues based on this idea were addressed in the literatur^. A 
formulation if IR quantum gravity in curved space was suggested [15] and the logarithmic 
corrections to scaling relations in the IR regime were studied [16]. The IR dynamics of the 
conformal factor was also investigated in four-dimensional quantum gravity with torsion [17j 
and a possible curvature induced phase transitions in IR quantum gravity was suggested 

m- 

In this paper we investigate the effective low energy gravity supplemented by the Bianchi 
identity constraint and its implications for cosmology. We obtain a model in which, besides 
a scale dependent cosmological constant, dark matter (DM) has a noncanonical cosmological 
scale dependence. We find that, under reasonable assumptions in the cosmological context, 
the running of the DM particle mass is phenomenologically acceptable. We demonstrate that 
a reasonable cosmology is obtained for a range of parameters required to fit the observed 
spectral index of the primordial density fluctuations. 

We organize the paper as follows. In section [2] we briefly discuss the large scale effects of 
the trace anomaly. In section [3] we introduce a cosmological model based on a generalized 
Bianchi identity. Finally, in section H] we summarize our results and give concluding remarks. 



2 Trace anomaly and effective low energy gravity 

In this section we summarize the basic ideas and results of Antoniadis, Mazur, and Mottola 
[T9] concerning the gravitational trace anomaly. The effects of the trace anomaly in the 
conformal sector of gravity may be studied by making use of the conformal parameterization 

9,uix) = e'-^^^g,,ix), (1) 

where is a fixed fiducial metric. The total low energy effective action is 

•S* = (SeH + 'S'niatt + 'S'anom; (2) 

where 

is the classical Einstein-Hilbert action, Smatt is the action that contains the matter fields and 
5'anom IS the auomaly induced effective action [131 1201 [211 122] 

(4) 



bFa + b' [E- ^Di? ) a + 2h'a^ia 



5'anom J ^ 

The differential operator A4 defined as 



A4 = + 2W^^^^, - hiu + hy^R)^^ , (5) 



^For recent reviews, see [T^ 114) . 
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is the unique conformally invariant 4-th order operator. The parameters b and b' are the 
coefficient that muhiply respectively the square of the Weyl tensor 



F = C.^.xC^"'^ = R.upxR^"'' - 2i?^.i?^^ + ^R' (6) 
and the Gauss-Bonnet term 

E = R^upxR^""'' - AR^^R^" + R^ (7) 

that appear in the most general expression for the 4-dimensional trace anomaly [23]. These 
parameters depend on the matter content of the theory coupled to cr. If only the contribution 
of free massles fields are taken into account, b and b' take on the values 

^ ^ ^—{Ns + ?>Nf + l2Nv - 8) + fegrav 
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where Ns, Np, and Ny are the numbers of scalar, Weyl fermion and vector fields, respectively. 
The integers —8 and —28 come from the spin-0 metric and ghost fields while fegrav and b'^^..^^ 
are the contributions of the spin-2 metric fields. In the following we treat as a free 
parameter since the contributions beyond the Standard Model are not well known and the 
spin-2 contributions in and depend on the model of quantum gravity and are still an 
open issue. The one loop calculations in the Einstein theory and in the Weyl-squared theory 
give similar results [2Tj. However, it is not obvious that the Weyl-squared Lagrangian is 
appropriate to account for the contribution of traceless graviton degrees of freedom and the 
Einstein theory is plagued with several well-known difficulties. Although it follows from (Q 
that > for all free matter fields, we allow to take on negative values. Negative sign 
contributions can be obtained, e.g., in some extended models of conformal supergravity [24j 
(for additional references, see also [Hj). 

The scale invariant effective theory that gives rise to (jl]) has a non-trivial IR dynamics 
owing to the existence of a non-trivial IR stable fixed point [LSI- The scale invariance in 
this theory persists even at the quantum level. The sectors of a theory, the scale invariance 
of which persists at the quantum level, have recently been dubbed "unparticle stuff' [25J. 
These sectors, if coupled to the standard model sector, seem to cause novel observable effects 
which could perhaps be detected in the future experiments at TeV scale. In particular, one 
could also expect that the unparticle stuff gives additional contributions to Q^. 

The quantum fluctuations of the conformal factor are responsible for a screening of the 
cosmological and inverse Newton coupling constants fTZl [TU] . The anomalous dimension of 
an operator with the classical conformal weight u is given by the quadratic equation [T^l 



/5p = 4-p+^, (10) 



with the solution 

f3, = Q'-^/Q'-i8-2p)Q\ (11) 



3 



where p = 4 — is the classical conformal codimension. The full scaling dimension Ap is 
related to the classical dimension by ^7] 

A.^4(l-|). (12) 

The operators 5*0 and 5*2 appearing in ([3]) aquire anomalous dimensions Pq and P2, respec- 
tively, and scale with volume according to 

^0 ~ V^; ^2 ~ V^'/^\ (13) 

whereas the corresponding couplings scale inversely 

^ y-l; 1 (14) 



SnG ' 167rG 

By similar considerations one finds the scaling laws for fermion and boson masses 

mp ~ -^-ftZ/Jo. ^1 _ y-MPo^ (15) 

Another effect of the quantum fluctuations of the conformal factor concerns the departure 
of the fractal space-time dimension from its classical value [19]. It turns out that the volume 
V does not scale with geodesic distance / naively as V l^. Rather, it scales according to 

V r^f^ , (16) 

where d-a is the Hausdorff dimension which classically equals 4. The calculation based on the 
quantum gravity distance defined by the heat kernel of the operator A4 yields the Hausdorff 
dimension expressed in terms of the parameter Q"^ 



du = (17) 

where jSs and Po are given by ( fTTll . In the course of the cosmological evolution the physical 
geodesic distance / scales as / ~ a and hence the volume scales with a as 

V ~ a'^^. (18) 

As has been emphasized in [T9j, the scaling relations f|T4l) and f|T5l) of dimensionful quantities 
are not directly physically relevant since the units in which the volume is measured have not 
been specified. A physically meaningful scaling relation is obtained when a product of 
powers of two quantities is formed so that its naive dimension is zero. In this way one of the 
quantities is measured in units of the other. 

Combining (IT^ . (fTSjl . and (fTSll . the net effect is a cosmological scale dependence of the 
dimensionless quantities [19] 

GA = GoAoa'^, (19) 
Gml = GomlQa", (20) 
Gml = Gomiy, (21) 
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where the exponents jj, and u are given by 



= c?H 2— - 1 ; z/ = ^ • (22) 

\ Po / Po 

In the hmit of large one finds 

Thus, for positive the cosmological constant decreases, the fermion masses increase, and 
the boson masses remain constant with increasing cosmological scale a, when these quantities 
are measured in units of the Planck mass. 



3 Cosmology by the generalized Bianchi identity 

A scale-setting procedure based on the implementation of the generalized Bianchi identity 
was established and successively applied J28j both to the effective field theory of gravity and 
matter, and to the nonperturbative quantum gravity [29]. Besides, it was found [30l[3T] that 
both theories are consistent with holographic dark energy [32], provided the running scale 
was identified with an infrared cutoff roughly equal to the inverse size of the system. 

It seems reasonable to assume that the quantum effects of the conformal factor which 
we have discussed in section |2] do not alter the Einstein field equations at large distances, 
apart from the gravitational dressing of G and A due to these effects. Then, the contracted 
Bianchi identity of the Einstein tensor yields the conservation law 

W[G{T,, + g^,PA)]=0, (24) 

where the energy- momentum tensor takes the usual perfect fluid form T^^, = {p + p)u^Uy + 
pg^u. In the comoving frame with FRW metric, equation becomes 

[G{p + pA)] + 3Gip + p) = 0. (25) 

da 



where the cosmological scale a satisfies the Friedmann equation in fiat space 

H'^i-] =^(p + Pa). (26) 



aj 3 

In fl25|) we implicitly assume scale dependent pA and G and the scale dependence of the DM 
density p need not be canonical a~^~^^. Furthermore, we assume that matter is nonrela- 
tivistic, i.e., that the matter energy density can be written as p = ^jnjmj, where rij are 
the particle number densities and rrii are the masses of the particle species. If nonrelativistic 
matter consists of / fermionic and k bosonic species, we have I + k + 2 equations (equations 
(fT9|) - (!2T]) and (p5|) ) for 2 + 21 + 2k quantities. Although we do not expect k and / to be much 
larger than 1, additional assumptions are needed in order to solve the equations uniquely. 
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Postulating (125|) . the scale behavior of p depends on the scaling of A and G. Equation 
( IT9|) gives the scaling of the dimesionless product AG and does not determine the scaling 
of the dimensionful parameters G and A separately. However, we are allowed to choose the 
units of measurement such that one chosen dimensionful parameter is fixed. We may, e.g., 
investigate three obvious choices when one of the three quantities G, A, or p\ is fixed. 

A more general case that includes the above mentioned choices is obtained if we alow G 
and A to vary as powers of a restricted only by equation f|T9|) . Hence, we set 



G = Goa"; A = Aoa^"", (27) 

where a is an arbitrary parameter. With this ansatz, p may be expressed in terms of the 
fixed dimensionful parameter A^/G'^"" = Aq/Gq"" 



P(«) = [q^) = ^APcfia), (28) 

where /(a) = 87r(GoAo)^"^'^~^p*(a) is a dimensionless function of a, pc is the critical density 
at present, and the constant Q\ which may be fixed from observations is of the order of the 
present fraction of DE density. 

Plugging fl27|) and fl28|l in fl25|) and neglecting the DM pressure we obtain a differential 
equation for / 

+ (3 + a)f + (p- a)a^~2a _ g, (29) 

da 

with the solution 

/ = Ga-=^-" a^-'". (30) 

3 + /i — a 

The integration constant G is for small p and a roughly G ~ (1 — Qa)/Qx so that p fits 
the present DM density. Although the parameter a is arbitrary, the small values -C 1 
and |a| -C 1 are phenomeno logically desirable since the variation of A and G should not 
be too large to spoil the observational constraints. Hence, equation (150]) implies a slight 
modification of the DM density. However, it may be easily seen that the effective DM density 
in the Friedmann equations remains canonical. Using fl27j) . fl28l) . and fl30|) the Friedmann 
equations may be written in the form 

2 SvrG 

(Pdm + Pde), (31) 



a / 3 

(Pdm + Pde + 3pde) , (32) 



d 47rGo 



a 3 

where the effective DM and DE densities are given by 

Pdm = (1 - flT>E)Pca~^, (33) 
PDE = f^DEPca^"", (34) 

and the effective DE pressure is 

PDE = -(l + ^] PDE . (35) 
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In (l33l) and we have introduced the constant Qbe which we identify with the present 
fraction of DE. In this way we fix the arbitrary constants Q/<^ and C which are related to 
Qbe by 

n» = (l + ^i^) !Joe; (36) 

Equations (l3Tl) -( l35l) show that the models satisfying (|T9l) and (1271) closely mimic the cosmol- 
ogy with standard cold DM and dark energy with a constant equation of state flM|) (XCDM 
cosmologies [33]), at least at the level of the global evolution of the universe. There is an 
obvious twofold implication of this result. First, in the analysis of the observational data, 
especially the SN la observations, one should bare in mind that a good fit to an XCDM 
cosmology may be interpreted as a signal for the dynamics given by fl27|) . Second, the ob- 
servational constraints on XCDM cosmologies can be used to estimate the parameters of 
the models defined by (fT9|) and (1271) . The values of a equal to 0, /x, and /i/2 correspond to 
fixed G, A, and pA, respectively. Although the simplest result is obtained for a = p, the 
natural choice is a = 0, i.e., fixed G, since a variation of G is subject to the most restrictive 
observational constraint^. In this case we obtain the DE equation of state 

^ = -l-|. (37) 
The observational constraint on the equation of state [3] w = —0.97 ± 0.07 yields 

-0.3 < ;u < 0.12, (38) 
which in turn constrains the parameter Q^. Using (l23ll we obtain 

< -33 or > 13 , (39) 

where < would imply a phantom cosmology. 

Another constraint on comes from the primordial density fluctuations [26]. It is usual 
to assume that the two point correlation function of the primordial density fluctuations 6p, 
should behave like 

{6p{k)6p{-k)) k"" , (40) 

where the exponent n, called the spectral index, need not be constant over the entire range of 
wave numbers. According to Harrison and Zel'dovich [35] the primordial density fluctuations 
should be characterized by a spectral index n = 1. In other words, the observable giving rise 
to these fluctuations has naive scaling dimension p = 2. This naive scaling dimension reflects 
the fact that the density fluctuations are related to the metric fluctuations by Einstein's 
equations SR ~ G6p in which the scalar curvature is second order in derivatives of the 
metric. Hence, the two point spatial correlations {6R{x)6R{y)) should behave as \x — y\~^ 
or in Fourier space. 

More generally, as a consequence of conformal invariance the two-point correlation func- 
tion of an observable O with dimension A is given by 

{0{x)0{y))r.\x-y\-'^, (41) 



^For an extensive discussion of models with variable A and fixed G, see 
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at equal times in three dimensional flat space. In fc-space this becomes 



{0{k)0{-k)) ~ \k\ 



2A-3 



(42) 



Thus, the spectral index of an observable of dimension A is defined by 



n = 2A - 3 



(43) 



Hence, the Harrison-Zel'dovich spectral index n = 1 corresponds to the classical dimension 
A = 2 of the primordial density fluctuation. 

If the conformal fixed point behavior [13] described in section [2] dominates at cosmological 
scales then the scaling dimension Ap of an observable with classical dimension p is given by 
( IT2l) as required by the conformally invariant fixed point for gravity. As a consequence 
of (H3l) . with A2 instead of A = 2, a deviation from the Harrison-Zel'dovich spectrum is 
obtained. For large one finds [121 EE] 



A comparison of n thus calculated with recent observations yields a constraint > 80. 

The favored WMAP value seems to be n = 0.95 which requires a large negative ~ —80. 
With this value, we obtain dark energy of the phantom type with w ~ —1.02, which is 
consistent with SN la and WMAP observations. This result justifies a relaxation of the 
allowed range for the parameter Q^. The above considerations show that using a single 
value of the parameter in our model it is possible to satisfy the observational constraints 
for two essentially unrelated phenomena: the present accelerated expansion of the universe 
( !37l) and the spectral index of primordial density fluctuations (l44l) . The required negative 
value of cannot be easily accommodated within the framework of the present theory, 
but the phenomenological potential of negative is a clear incentive to search for new 
mechanisms which could bring into the negative realm. 

The ansatz fl27|) is not the most general and does not cover all interesting possibilities. 
For example, it does not include a natural starting assumption that the total energy density 
of nonrelativistic matter scales canonically with a. With this assumption, equations flTIJl) 
and fl2^ fully determine the evolution of A and G. However, the canonical scaling of the 
matter energy density, p ~ a~^, combined with the scalings (l20l) and (!2T!) . implies that 
the particle number densities Ui no longer vary as and hence, the particle numbers of 
individual species are not conserved. Although the assumption of canonical scaling of the 
matter density yields a mathematically consistent model, we find a strong disagreement with 
observations in a wide parametric range. In particular, for negative fi we obtain a maximal 
allowed value of the redshift of the order of 1 which is clearly not acceptable. 

Another interesting model, not covered by fl27|) . is obtained if one assumes that all relevant 
particle species are fermions and that the corresponding particle number densities scale 
canonically, i.e., rij ~ a^'^. In this case one can express the mass in terms of G and a from 
fl2U]) . A in terms of G and a from flTI?]) and solving fl2^ obtain an evolution equation for G 
as a function of a. However, numerical solutions of this equation show that the resulting 
scaling of G with a is too strong to satisfy observational bounds, even for small values of 
parameters fj, and u. 




(44) 



8 



4 Conclusion 



We have studied some DE and DM aspects of the low energy effective theory of gravity. This 
theory is a modified general relativity in which the Einstein-Hilbert action is supplemented 
with the non-local terms induced by the trace anomaly of massless fields. These nonlocal 
terms do not decouple for scales E -C Mpi and therefore become increasingly important for 
the present and future time cosmology. The part of the action that stems from the trace 
anomaly contains all infrared relevant terms which are not contained in the local action. 

The testing of the theory vs. observational astrophysics and cosmology is a long term 
project, which includes the use of as a dynamical source for Einstein's equations. 

However, we believe that the effective theory with running G and A supplemented with the 
generalized Bianchi identity may be successfully confronted with cosmological observations. 

The effective low energy gravity in the conformal sector predicts a cosmological scale 
dependence of various dimensionless quantities (equations ( TT9|) -(l2Tl)). The scale dependence 
of DM and DE densities, being dimensionful quantities, depends on the choice of a fixed 
dimensionful parameter. In particular, fixing the Newton constant G yields a cosmological 
constant scaling with a as A ~ a'^ and a noncanonical scaling of the DM energy density given 
by (l30l) with a = 0. The effective DM and DE densites yield a reasonable cosmology if the 
exponent fi is small and restricted by the constraint (155]) . This constraint in turn yields a 
constraint (1391) on consistent with the observational bounds on the spectral index of the 
primordial density fluctuations. 

In our approach the parameter that appears in the effective action (jl]) induced by the 
gravitational anomaly is treated as a free parameter. This parameter could, in principle, be 
calculated if one had a complete information on the conformally invariant sector. Unfortu- 
nately this sector is yet not well known so a precise theoretical prediction for remains an 
open problem. 

Acknowledgments 

This work was supported by the Ministry of Science, Education and Sport of the Repub- 
lic of Croatia under contracts No. 098-0982930-2864 (N.B., E.G., H.N., H.S.) and No. 
098-0982887-2872 (R.H.), and partially supported through the Agreement between the As- 
trophysical Sector, S.I.S.S.A., and the Particle Physics and Cosmology Group, RBI. 

References 

[1] A.G. Reiss et al., Astronom. J. 116 (1998) 1009; 
S. Perlmutter et al., Astrophys. J. 517 (1999) 565; 

S. Perlmutter, M.S. T urner, M.J. White, Phys. Rev. Lett. 83 (1999) 670 
|arXiv:astro-p h/990l052l . 

[2] C.L. Bennett et al., Astrophys. J. Suppl. 148 (2003) 1 |arXiv:astro-ph/0302207l; 
D.N. Spergel et al., Astrophys. J. Suppl. 148 (2003) 175 |arXiv:astro-ph/0302209]; 
H.V.P. Peiris et al, Astrophys. J. Suppl. 148 (2003) 213 [arXiv:astro-ph/0302225] . 



9 



[3] D.N. Spergel et al., Astrophys. J. Suppl. 170 (2007) 377 |arXiv:astro-ph/0603449| . 

[4] T. Padmanabhan, Phys. Kept. 380 (2003) 235 farXiv:he p-th/0212290| ; 
P.J.E. Peebles, B. Ratra, Rev. Mod. Phys. 75 (2003) 559; 
S. Nobbenhuis, Ph D. Thesis, arXiv:gr-qc/0609011[ and references therein. 

[5] S. Weinberg, Physica (Amsterdam) 96A (1979) 327. 



[6] J.F. Donoghue, arXiv:gr-qc/9712070; 

J.F. Donoghue, Helv. Physica Acta 69 (1996) 269 farXiv:gr-qc/9607039|; 
J.F. Donoghue, Phys. Rev. D 50 (1994) 3874 |arXiv:gr-qc/9405 057|. 

[7] N.D. Birell, P.C.W. Davies, Quantum Fields in Curved Space, Cambridge University 
Press, Cambridge, 1992. 

I.L. Buchbinder, S.D. Odintsov, I.L. Shapiro, Effective Action in Quantum Gravity, 
lOP, Bristol, 1992. 

[9] I. Shapiro, J. Sola, Phys. Lett. B 475 (2000) 236; 
I. Shapiro, J. Sola, JHEP 0202 (2002) 006. 

[10] A. Babic, B. Guberina, R. Horvat, H. Stefancic, Phys. Rev. D 65 (2002) 085002. 

[11] I.L. Shapiro, J. Sola, C. Espana-Bonet, P. Ruiz-Lapuente, Phys. Lett. B 574 (2003) 149; 
F. Bauer, Class. Quant. Grav. 22 (2005) 3533; 
LL. Shapiro, J. Sola, H. Stefancic, JCAP 0501 (2005) 012; 
J. Sola, H. Stefancic, Mod. Phys. Lett. A 21 (2006) 479; 
J. Grande, J. Sola, H. Stefancic, JCAP 0608 (2006) Oil; 
J.C. Fabris, I.L. Shapiro, J. Sola, JCAP 0702 (2007) 016. 

[12] I. Antoniadis, P.O. Mazur, E. Mottola, New J. Phys. 9 (2007) 11 |arXiv:gr-qc/0612068| . 

[13] I. Antoniadis, E. Mottola, Phys. Rev. D 45 (1992) 2013. 

[14] I.L. Shapiro, |arXiv:hep-th/061016 8l 

[15] S.D. Odintsov, Z. Phys. C 54 (1992) 531. 

[16] I. Antoniadis, S.D. Odintsov, Phys. Lett. B 343 (1995) 76 |arXiv:hep-th/9411012| . 

[17] I. Antoniadis, S.D. Odintsov, Mod. Phys. Lett. A 8 (1993) 97 9 |arXiv:hep-th/92 12140|; 
E. Elizalde, S.D. Odintsov, Int. J. Mod. Phys. D 2 (1993) 51 |arXiv:hep-th7920 5048 ] . 

[18] E. Elizalde, S.D. Odintsov, Mod. Phys. Lett. A 8 (1993) 3325 |arXiv:hep-th79211008] ; 
Phys. Lett. B 315 (1993) 245 | arXiv:hep-th/93 08019|. 

[19] I. Antoniadis, P.O. Mazur, E. Mottola, Phys. Lett. B 444 (1998) 284 
I arXiv:hep-th/9808070 | . 

[20] R.J. Riegert, Phys. Lett. B 134 (1984) 56; 

E.S. Fradkin, A.A. Tseytlin, Phys. Lett. B 134 (1984) 187. 



10 



[21] I. Antoniad is, P. O. Mazur, E. Mottola, Nucl. Phys. B 388 (1992) 627 
|arXiv:hiH ?h/9205015|. 

[22] I . Antoniadis, P.O. Mazur, E. Mottola, Phys. Rev. D 55 (1997) 4770 
|| arXiv:hep-th79509 169| . 

[23] M.J. Duff, Nucl. Phys. B 125 (1977) 334. 

[24] E.S. Fradkin, A.A. Tseytlin, Phys. Reports 119 (1985) 233. 

[25] H. Georgi, Phys. Rev. Lett. 98 (2007) 221601 (arXiv:hep-ph/0703260). 

[26] 1. Antoniadis, P.O. Mazur, E. Mottola, Phys. Rev. Lett. 79 (1997) 14 
| arXiv:astro-ph/9611208] . 

[27] I. Antoniadis, P.O. Mazur, E. Mottola, Phys. Lett. B 394 (1997) 49 
| arXiv:hep-th/9611145j . 

[28] A. Bab ic, B. Guberina, R. Horvat, H. Stefancic, Phys. Rev. D 71 (2005) 124041 
^arXiv: astro-ph /040 75 72] . 

[29] M. Renter, Phys. Rev. D 57 (1998) 971; 

A. Bonnano, M. Renter, Phys. Rev. D 65 (2002) 043508; 
A. Bonnano, M. Renter, Phys. Lett. B 527 (2002) 9. 



[30] R. Horvat, Phys. Rev. D 70 (2004) 087301 [arXiv :astro- ph/0404204] . 
[31] B. Guberina, R. Horvat, H. Stefancic, JCAP 05 (2005) 001. 

[32] A. Cohen, D. Kaplan, A. Nelson, Phys. Rev. Lett. 82 (1999) 4971 
[arXiv:hep-th/9803132|. 



[33] M.S. Turner, M.J. White, Phys. Rev. D 56 (1997) 4439 | arXiv:astro-ph/9701138| . 

[34] J.M. Overduin, F.I. Cooperstock, Phys. Rev. D 58 (1998) 043506. 

[35] E.R. Harrison, Phys. Rev. D 1 (1970) 2726; 
Ya.B. Zel'sovich, MNRAS 160 (1972) IP 



11 



